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ABSTRACT 

We report the results of monitoring of V4334 Sgr (Sakurai's Object) at 450 fim and 
850 /xm with SCUBA on the James Clerk Maxwell Telescope. The flux density at both 
wavelengths has increased dramatically since 2001, and is consistent with continued 
cooling of the dust shell in which Sakurai's Object is still enshrouded, and which still 
dominates the near-infrared emission. Assuming that the dust shell is optically thin at 
sub-millimetre wavelengths and optically thick in the near-infrared, the sub-millimetre 
data imply a mass-loss rate during 2003 of '^3.4 ± 0.2 x 10~^ Mq yr~^ for a gas-to- 
dust ratio of 75. This is consistent with the evidence from 1 — 5 /im observations that 
the mass-loss is steadily increasing. 
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1 INTRODUCTION 

V4334 Sgr (Sakurai's Object) continues to be the subject of 
intense interest. It is a low- mass (^O.TM©) star that is re- 
tracing its evolution along the Hertzsprung-Russell diagram 
following a very late thermal pulse (Herwig 2001). The star 
is at the centre of a faint planetary nebula (PN) of diam- 
eter 40" (e.g. PoUacco 1999; Bond & PoUacco 2002) and, 
after an early incarnation as an F-type star, it underwent a 
transformation to a carbon star in 1997. 

Little is known observationally about the evolution of 
the star since mid- 1998, when it was completely obscured by 
an optically thick dust shell from which (as of April 2004) it 
has still not emerged; a dramatic illustration of the obscu- 
ration is given by Bond & PoUacco (2002). Since the major 
dust event of 1998 the 1 — 5 /im spectum is black-body-like 
and is simply that of the dust shell. Analysis of the post- 
1998 infrared (IR) observations has shown that the dust is 
carbon, primarily in amorphous form (Geballe et al. 2002; 
Tyne et al. 2002; Eyres et al., in preparation). The 1 — 5 /xm 
spectra during the period 1999-2001 were consistent with 
a mass-loss that increased from ~ 1.8 x 10~^ M0 yr~^ to 
~ 4.4 X 1O~®M0 yr~^ (for the gas-to-dust ratio of 75 we 
adopt here; see §3 below), while the maximum grain radius 
in a grain size distribution n{a) da oc a^' da {q ~ 3) in- 
creased by a factor ~ 3 (Tyne et al. 2002). 

We first observed V4334 Sgr at sub-millimetre wave- 



lengths with the SCUBA instrument (Holland et al. 1999) on 
the 15-m James Clerk Maxwell Telescope (JCMT) in 2001 
August, and reported marginal detections at both 450 /im 
and 850 /im (Evans et al. 2002, hereafter Paper I). 

We report here the results of our continued monitoring 
of Sakurai's Object at these wavelengths. 



2 OBSERVATIONS AND DATA REDUCTION 

We have observed Sakurai's Object on the dates given in 
Table 1, in which we include our 2001 August observation 
from Paper I. On each occasion atmospheric extinction was 
regularly monitored by taking skydip measurements, and 
pointing sources in the observatory catalogue were used to 
check telescope pointing and focussing; Uranus and Mars 
were used as calibrators. The beamwidths at 450 fim and 
850 /im are 7" and 14" respectively. 

SCUBA was used in photometry mode, which provides 
simultaneous observations at 450 /im and 850 /im; the obser- 
vations were carried out in good (sub-millimetre) weather 
conditions to ensure that data at both 450 /im and 850 /im 
were obtained. However, as is well-known (Holland et al. 
1999), radiation at the shorter wavelength is severely ef- 
fected by atmospheric water vapour, and hence detection 
and calibration is always more problematic at 450 /im even 
when weather conditions are favourable. The observations 
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Table 1. Sub-millimetre flux densities fu of V4334 Sgr. Errors are la. 



Date 
D/M/Y 


MJD 
- 2450000 


FCF (m 
450 /im 


Jy mV-i) 
850 Atm 


(mjy) 


Detection 
level (it) 


(mJy) 


Detection 
level (a) 


Comment 


18/08/2001 


2140 


309 ± 72 


219 ± 30 


96.6 ± 37.5 


3.2 


10.6 ±3.4 


3.5 


Paper I 


09/01/2003 


2649 


378 ± 49 


193 ± 5 


55.2 ± 14.9 


4.2 


35.7 ± 2.4 


17 




05/05/2003 


2765 


322 ± 3 


202 ± 1 


157.6 ±25.8 


4.1 


52.1 ± 2.7 


20 




15/06/2003 


2806 


440 ± 5 


217 ± 2 


213.2 ±29.1 


7.4 


61.3 ±2.9 


28 




16/06/2003 


2807 


401 ± 2 


213 ± 1 


113 ±41 


2.8 


58.3 ±3.2 


18 
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Figure 1. The 450 nm (left) and 850 /im (right) light curves of Sakurai's Object. 



consisted of 2 (or 4) sets of 50 integrations, to give total 
on-source integration time of 30 (or 60) minutes. First order 
sky subtraction was achieved by chopping 60" in azimuth. 

The data were reduced using the SURF package (Jen- 
ness & Lightfoot 2000). The frames were flat-fielded and 
corrected for atmospheric extinction, and the outer bolome- 
ters in the array were used for second order sky subtraction; 
'spikes' in the data were clipped to remove any points fur- 
ther than 3a from the mean. The flux conversion factors at 
850 fjxn and 450 /itm, together with the detection levels (the 
signal-to-noise ratios in the measured signal) for Sakurai's 
Object are listed in Table 1. 

Sakurai's Object was strongly detected at both 450 /im 
and 850 /xm, and the flux densities are listed in Table 1. 
As usual, a significant fraction of the uncertainty in the 
flux densities for Sakurai's Object listed in Table 1 arises 
from the uncertainty in the flux calibration, particularly at 
450 /im. In the Tabic the flux errors include both statisti- 
cal errors and calibration uncertainties. The measurement 
at 450 fj,m for 2003 June 16 is marginal and is not used in 
what follows. 



3 DISCUSSION 

The 450 /im and 850 /im light curves are shown in Fig. 1. 
We see a consistent and significant rise in the 850 /im flux 
density between 2001 August and 2003 Sepember; the be- 
haviour of the 450 yum data is consistent with the rise at 
850 (Um, but of course the 450 /xm calibration is less certain. 
In Paper I we considered the possibility that the 450/840 fim 
emission detected in 2001 August might have originated in 
the old PN but we concluded that the emission originated 
in material ejected in the 1995 event. The fact that the sub- 



millimetre fiux densities rise so dramatically is consistent 
with this conclusion. 

We first note that the sub-millimetre data may imply a 
drop in the /3- index of the dust (defined in the usual way such 
that the dust emissivity at frequency u is <x v'^) , from ~ 1.5 
in 2001 August to ~ in 2003. This may indicate changes 
in (a) the nature of the dust (e.g. as grains grow, or are 
processed) and/or (b) the dust distribution, allowing hotter 
grains (which have fiatter emissivities, cf. Mennella et al. 
1998) to contribute more to the observed fiux. However any 
conclusion along these lines must await detailed modelling 
combining both UKIRT (e.g. Tyne et al. 2002, Eyres et al., 
in preparation) and JCMT datasets; such an analysis will 
be reported elsewhere. 

The referee has pointed out that the sub-millimetre rise 
could be due to the anticipated rise in effective temperature 
as the star retraces its post-AGB evolution towards the PN 
phase, and the consequent heating of the inner dust shell. 
While observational evidence that this process may have 
started has been presented by Kerber et al. (2002), we do not 
consider that this accounts for the rise in the sub- millimetre 
fiux for reasons which may be discerned from Fig. 2. 

In this Figure we plot the JCMT data from 15 June 
2003 (Table 1), together with the 2-5 /xm data obtained in 
September 2003 from our UKIRT programme (Eyres et al., 
in preparation); we assume that variations are sufficiently 
slow that the data can be combined. Both sets of data are 
reasonably fitted by a ~ 360 ± 10 K black-body (cf. Equa- 
tion (2) below, and 350 ± 30 K for the UKIRT data alone; 
Eyres et al. 2004), suggesting that the 2-5 /nm and sub- 
millimetre emission have common origin. 

We might envisage a scenario in which the sub- 
millimetre emission in Fig. 2 comes from inner, hotter, dust 
whose near-IR emission does not penetrate the outer region 
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Table 2. Properties of the dust in the environment of Sakurai's Object. A gas-to-dust ratio of 75 is assumed; see text for details. 
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Date 
D/M/Y 


MJD 
-2450000 


Md Td (10- 


-5Mq K) 


Td (K) 


Md (10- 


'M0) 


Me, (10- 




450 /im 


850 Aim 


450 ^lm 


850 Aim 


450 /j.m 


850 Mm 


18/08/2001 


2140 


4.6 ±1.8 


3.3 ± 1.1 


504 


0.9 ±0.4 


0.7 ±0.2 


0.7 


0.5 


09/01/2003 


2649 


2.7 ±0.7 


11.1±0.8 


418 


0.6 ±0.2 


2.7 ±0.3 


0.5 


2.0 


05/05/2003 


2765 


7.5 ± 1.2 


16.2 ±0.8 


401 


1.9 ±0.4 


4.0 ±0.5 


1.4 


3.0 


15/06/2003 


2806 


10.2 ± 1.4 


19.1 ±0.9 


395 


2.6 ±0.4 


4.8 ±0.5 


1.9 


3.6 


16/06/2003 


2807 




18.2 ± 1.0 


395 




4.6 ±0.5 




3.5 




A [/im] 

Figure 2. Combined UKIRT/JCMT data for Sakurai's Object in 
late 2003 to give the 2.5-850 fim spectral energy distribution; the 
IJKIRT data have been binned for clarity. The curve is a 360 K 
black-body, fitted by eye. See text for details. 
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Figure 3. Increase in circumstcUar mass around Sakurai's Object 
as deduced from 450 fim data (squares) and 850 /im data (trian- 
gles), for a gas-to-dust ratio of 75. Lines are included to guide the 
eye. See text for details 



of the dust shell, while the observed ncar-IR emission comes 
from dust in the outer shell having very steep /3-index. How- 
ever such a situation would be ad hoc and contrived. Indeed 
that the /?-index is close to zero in late 2003 is supported 
both by the 450/850 flux ratio, and by the closeness of the 
2-850 /Lim data to a black-body at ~ 360 K (see Fig. 2). 

We therefore make the assumptions that (i) the emis- 
sion we detect at 450 /xm and 850 iim is from the dust shell 
(any emission at these wavelengths from the star - what- 
ever its effective temperature - is at the few /iJy level at 
most), (ii) the dust shell is optically thin at sub- millimetre 
wavelengths and (iii) the sub-millimetre emission is on the 
Ray leigh- Jeans tail. 

Wc justify (ii) by noting that the optical depth r at 
wavelength A (in fiin) is ~ 0.55 t(0.55 /im)/A(/im) <g 1, 
since the optical depth in the visual is ~ 10 (see Tyne 
et al. 2002), and (iii) is valid if the dust temperature is 
> 200 - 300 K. Eyres ct al. (2004) find that a 1 - 5 /im spec- 
trum of Sakurai's Object obtained on 2003 September 8 is 
well-described by a black-body at temperature ~ 350±30 K, 
consistent with continued cooling of the dust shell; thus 
(iii) is certainly justified. We therefore have that 
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where the constant A = 9.60 x 10"^ for 450 /im, and 6.23 x 
10~® for 850 /tm. In Equation (1) Mj is the dust mass, D is 
the distance of Sakuai's Object (see Kimeswenger 2002), Td 



is an appropriate value for the dust temperature (see below) , 
K is the absorption coefficient for the dust, and 7 depends 
on physical conditions in the dust shell (e.g. 7 s 1 for an 
optically thin, isothermal shell). 

Although Equation (1) applies to an optically thin, 
isothermal shell, it is also accurate, to within a factor 
~ 2/(1 + {Ri/RoY^'^) < 2 (where Ri,Ro are respectively 
the inner and outer radii of the dust shell), to the case of 
a dust shell with grain density law n(r) cx that is opti- 
cally thick at short (~ optical/ncar-IR) wavelengths, opti- 
cally thin in the sub-millimetre, and for which the dust tem- 
perature depends on distance r from the star as Td oc r~^^'^ 
(see Ivezic & Elitzur 1997). In this case Td in Equation (1) 
is the temperature of the coolest dust; we take this to be 
the dust seen in the 1 — 5 /im data, which still detect the 
dust 'pseudophotosphere' (Tyne et al. 2002, Eyres et al., in 
preparation) . It is not unreasonable to assume that the dust 
shell around Sakurai's Object is geometrically thick, so that 
Ro/Ri ^ 1; for this case 7 = 2 and wc use this value here 
(sec also below). 

For the distance, we take D = 2 kpc (see Kimeswenger 
2002), and we have taken k from Mennella et al. (1998; their 
K values for amorphous carbon at 295 K, which is likely to 
be close to that appropriate for Sakurai's Object). 

We give values of Md Td, computed from Equation (1), 
in Table 2 for both wavelengths. For each epoch, we have 
estimated the temperature Td of the coolest dust by (i) fit- 
ting a black-body to 1 — 5 /tm data obtained in the course of 
our IR monitoring programme (see Tyne et al. 2002, Eyres 
et al. 2004, and paper in preparation) and (ii) fitting an 
exponential decline to the time-dependence of Td. We find: 
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Td (K) ~ 817 exp 



(815.5 - t') 



2735.3 



(2) 



where t' = MJD — 2450000 and the uncertainty in Td is typ- 
ically ^ 10% in reproducing the values in Tyne et al. (2002) 
and Eyres et al. (2004). We stress that this expression is 
not physically meaningful, and wc will refine it as wc obtain 
further IR data; we use it simply as a means of interpolat- 
ing the value of Td at the time of our scuba observations. 
These values arc included in Tabic 2, together with the cor- 
responding dust masses Aid] the uncertainties in Afd include 
an uncertainty of 10% in estimating Td. 

Neither the values of Md Td, nor of Md, should of course 
depend on wavelength, and the discrepancies between the 
450 /im and 850 /im values arc due to the uncertainties in 
the 450 /im calibration, the error in the interpolated tem- 
perature, uncertainty in the properties of the dust, and the 
fact that the 1 — 5/Ltm emission differs somewhat from a 
black-body. Nevertheless the data at both wavelengths are 
consistent with a steady increase in the dust mass between 
2001 August and 2003 June. 

The dust mass may be converted to a circumstellar mass 
McB if we assume a gas-to-dust ratio, and that the gas-to- 
dust ratio remains constant in time. The gas-to-dust ratio 
for solar abundances and in the interstellar medium is ~ 120 
(e.g. Whittet 2003), while a value 200 was assumed by Tyne 
et al. (2002). These values axe unlikely to be appropriate for 
the H-deficient, He- and C-rich wind of Sakurai's Object. If 
we take the abundances for 1996 July 3 from Asplund et al. 
(1999), assume that all O is locked up in CO (see Eyres ct 
al. 2004) and that the remaining C is tied up in dust, we get 
a gas-to-dust ratio of 26; this is most likely a lower limit. 

We assume an intermediate value of 75 here; the uncer- 
tainties in the gas-to-dust ratio and in the value of 7 clearly 
lead to a corresponding uncertainty in dust mass, but we 
believe that changes in the latter should be reliable. The 
corresponding circumstellar masses are given in Table 2; the 
dependence of circumstellar mass on time is shown in Fig. 3. 

Wc use these data to determine the average 
mass-loss rate over the period 2001 August to 2003 
June. Using the 850 /xm data, we find that M ~ 
1.6 ± 0.3 X 10~® Mq yr~^ over this period; the correspond- 
ing value from the 450 fim data is 4.9 ± 4.4 x 10"® Mq yr""^. 
Taking the 2003 observations only, the 450 /im and 850 /im 
data give M ~ 3.3 ± 0.3 x 10"^ M© yr'^ and M ~ 
3.6 ± 0.3 x 10~^ M© yr~^ respectively, with a weighted 
mean of M ~ 3.4 ± 0.2 x 10"^ M© yr'^ (cf. Fig. 3); for com- 
parison, the maximum mass-loss rate for post-AGB stars is 
- 1 X 10"" M© yr-i (van Winckel 2003). 

Tyne et al. (2002) used the dusty code (Ivezic & Elitzur 
1997) to fit the 1 — 5 /xm data with a carbon dust shell 
and, adjusted for a gas-to-dust ratio of 75, found that M 
increased from 1.8 x 10"® M© yr"^ to 4.4 x 10"® M© yr"^ 
over the period 1999 May to 2001 September. Our results 
suggest that the mass-loss from Sakurai's Object continues 
to rise. 



4 CONCLUSIONS 

We have presented 450 fjjn and 850 /ttm photometry of Saku- 
rai's Object, and find that the flux densities at both wave- 



lengths have increased steadily over the period 2001 Au- 
gust to 2003 June. These data, combined with dust temper- 
atures obtained from out IR monitoring programme, lead 
to a mass-loss in 2003 of ~ 3.4 ± 0.2 x 10"^ M© yr"\ con- 
sistent with a continuing rise since our IR observations of 
1999-2003. 

Monitoring of this object at sub-millimetre wavelengths 
is continuing and further observations will be presented else- 
where. 



ACKNOWLEDGMENTS 

Wc thank the referee, Stefan Kimcswcnger, for comments 
that helped to clarify some of the discussion. TRG is sup- 
ported by the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., 
on behalf of the international Gemini partnership of Ar- 
gentina, Australia, Brazil, Canada, Chile, the United King- 
dom, and the United States of America. The JCMT is op- 
erated by the Joint Astronomy Centre on behalf of the Par- 
ticle Physics and Astroonomy Research Council (PPARC) 
of the United Kingdom, the Netherlands Organization for 
Scientific Research and the National Research Council of 
Canada. We thank the PPARC Panel for the Allocation of 
Telescope Time for supporting of this project. Data reduc- 
tion was carried out using hardware and software provided 
by PPARC. 



REFERENCES 

Asplund M., Lambert D. L., Kipper T., Pollacco D., Shetrone M. 

D., 1999, A&A, 343, 507 
Bond H. E., Pollacco D., 2002, Keele Workshop on Sakurai's Ob- 
ject, Eds A. Evans & B. Smalley, Ap&SpSci, 279, 31 
Evans A., Geballe T. R., Tyne V. H., Pollacco D., Eyres S. P. S., 

Smalley B., 2002, MNRAS, 332, L69 (Paper I) 
Eyres S. P. S., Smalley B., Geballe T. R., Evans A., Asplund M., 

Tyne V.H., 1999, MNRAS, 307, Lll 
Eyres S. P. S., Geballe T. R., Tyne V. H., Evans A., Smalley B., 

Worters H. L., 2004, MNRAS, 350, L9 
Geballe T. R., Evans A., Smalley B., Tyne V. H., Eyres S. P. S., 

2002, Keele Workshop on Sakurai's Object, Eds A. Evans & 

B. Smalley, Ap&SpSci, 279, 39 
Herwig P., 2001, ApJ, 554, L71 
Holland W. S., et al., 1999, MNRAS, 303, 659 
Ivezic Z., Elitzur M., 1997, MNRAS, 8287, 799 
Jemiess T., Lightfoot J. P., 2000, Starlink User Note 216, Starlink 

Project, CLRC 

Kerber P., Pirzcal N., de Marco O., Asplund M., Clayton G. C., 

Rosa M. R., 2002, ApJ, 581, L39 
Kimeswenger S., 2002, Keele Workshop on Sakurai's Object, Eds 

A. Evans & B. Smalley, Ap&SpSci, 279, 79 
Mennella V., Brucato J. R., Colangeli L., Palumbo P., Rotundi 

A., Bussolctto E., 1998, ApJ, 496, 1058 
Pollacco D., 1999, MNRAS, 304, 127 

Tyne V. H., Evans A., Geballe T. R., Eyres S. P. S., Smalley B., 

Durbeck H. W., 2002, 334, 875 
van Winckel H., 2003, ARAA, 41, 391 

Whittet D. C. B.,2003, Dust in the Galactic environment, 2nd 
edition. Institute of Physics Publishing, Bristol & Philadel- 
phia 



Sub-mm evolution of Sakurai's Object 5 



This paper has been typeset from a KTgX file prepared 
by the author. 



